The Shift of the Baryon Acoustic Oscillation Scale: A Simple Physical Picture 
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A shift of the baryon acoustic oscillation (BAO) scale to smaller values than predicted by linear 
theory was observed in simulations. In this paper, we try to provide an intuitive physical under- 
standing of why this shift occurs, explaining in more pedagogical detail earlier perturbation theory 
calculations. We find that the shift is mainly due to the following physical effect. A measurement of 
the BAO scale is more sensitive to regions with long wavelength overdensities than underdensities, 
because (due to non-linear growth and bias) these overdense regions contain larger fluctuations and 
more tracers and hence contribute more to the total correlation function. In overdense regions the 
BAO scale shrinks because such regions locally behave as positively curved closed universes, and 
hence a smaller scale than predicted by linear theory is measured in the total correlation function. 
Other effects which also contribute to the shift are briefly discussed. We provide approximate ana- 
lytic expressions for the non-linear shift including a brief discussion of biased tracers, and note that 
the shifts are different in real and Fourier space due to a change of the shape of the BAO feature. 
We explain why reconstruction should entirely reverse the shift. Our expressions and flndings are 
in agreement with simulation results, and conflrm that non-linear shifts should not be problematic 
for next-generation BAO measurements. 
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I. INTRODUCTION 

The primordial photon-baryon plasma supports the 
propagation of sound waves, sourced by the initial den- 
sity fluctuations, until the photons and baryons decou- 
ple. At decoupling, a characteristic correlation length 
is hence imprinted into both the baryons and the pho- 
tons: the sound horizon at recombination. This man- 
ifests itself in the large scale structure of the universe 
as an excess probability of finding pairs of galaxies (or 
other tracers) at a separation equal to the sound horizon 
scale, known in this context as the Baryon Acoustic Os- 
cillation (BAO) scale. As the value of the sound horizon 
is well known from microwave background observations, 
the BAO scale (which is ~ 150 comoving Mpc) can be 
used as a "standard ruler" to infer the expansion history 
of the universe. It has been detected (either as a peak in 
the real space correlation function ^(r) or as a series of 
peaks in its Fourier transform, the power spectrum P{k)) 
by a number of large scale structure surveys (ll-[l5| , but 
its potential as a cosmological probe is nowhere near ex- 
hausted: due to its robustness to systematic errors, many 
future large scale structure surveys propose to measure 
the BAO scale to ~ 0.1% in order to infer the expansion 
history of the universe and the properties of dark energy 
to high precision, e.g. [l^. This proposed future research 
program relies on the BAO scale being a good standard 
ruler to the sub-percent level. 

However, the extent to which this is true needs to be 
examined carefully, as the BAO scale is somewhat mod- 
ified by non-linear gravitational efi^ects. Simulations of 
both dark matter and biased tracers [17, 23] show that 
the BAO standard ruler shrinks by - 0.3% [D{z)/D{0)f 
at redshift z, where D{z) is the growth factor. Ap- 
proaches usin g pe rturbation theory for both matter and 
biased tracers |18l - [2q | also find a shift of this order, which 
is explained as being due to the generation of out-of- 



phase oscillations due to non-linear mode coupling ef- 
fects. In this paper, we try to develop a more intu- 
itive physical picture for why the BAO scale shrinks due 
to nonlinearities. We first examine the effect of a long 
wavelength mode on the small scale correlation function 
and heuristically derive an approximate expression for 
the shift in the BAO scale. We then present a less ap- 
proximate pedagogical calculation of the shift with per- 
turbation theory. In the final sections, we discuss the 
effect of biased tracers and reconstruction on the shift of 
the BAO scale. 



II. THE PHYSICAL ORIGIN OF THE BAO 
SCALE SHIFT 

In this section, we will describe heuristically how the 
correlation function on small scales ^si^) (which contains 
a BAO-like feature) is affected by the presence of a mode 
6l with a wavelength larger than the scales of interest 
which are of order the BAO scale. We will obtain through 
this simple argument an expression for the shift in the 
BAO scale. In the following section, we will confirm and 
refine this expression using perturbation theory. 

As we are only interested in the correlation function on 
scales of order the BAO scale (and not much larger), 
we split the overdensity S into long wavelength perturba- 
tions 6l and short wavelength perturbations 6s so that 
S ~ Ss + Sl- The small scale correlation function can 
hence be defined as 



Ur) = {Ss{^ + v)Ss{^)), 



(1) 



We assume that the small scale correlation function 
is determined with a fixed background overdensity Sl 
in a small volume (smaller than the wavelength of the 
long mode, but larger than the BAO scale). As shown 
in f2l| and as can be deduced from Birkhoff's theorem, 
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physics in such a small volume is indistinguishable (to 
first order in S^) from physics in a universe with the same 
curvature as that induced locally by the overdensity Sl 
(and with the same early time limiting behavior as our 
flat universe). Hence the behavior of small scale density 
perturbations in the presence oi 6l is identical to their 
behavior in an appropriately curved universe. 

What are the effects of ^i,, and hence background cur- 
vature, on the small scale correlation function ^5? When 
Sl is positive in a volume, the volume behaves locally 
as slightly closed universe with positive curvature. The 
first effect of this is that all the features in the corre- 
lation function are contracted to smaller scales than in 
a flat universe by a factor F = acurved(i5L)/ao, where 
a is the scale factor. As the overall long wavelength 
density in this region must be simply related to F by 
1 + 5l = (j^)'^, the linear theory correlation function 
^so{r) becomes eso(r/F((5L)) « Cso((l + ^)f) due to this 
non- linear contraction (to first order in Sl, neglecting a 
second order term). 

The second effect of a long wavelength overdensity is 
on the growth of small scale perturbations. Short modes 
in regions where Sl is positive experience more growth 
than modes in regions where Sl is negative. This can 
again be explained as overdense regions locally behaving 
like a closed universe, in which the growth of local per- 
turbations (i.e. the growth of perturbations defined with 
respect to the background density of the closed universe 
pi = P + Spl) is enhanced by a factor 1 -I- 13/21(5l (see 
e.g. [2lj) and hence the local correlation function is en- 
hanced by(l + 13/21(5L)^- To convert a local correlation 
function defined with respect to pi to one defined glob- 
ally with respect to p, we must multiply it by (1 -f (5^)^. 
Hence the growth factor of the short wavelength corre- 
lation function due to the presence of a long wavelength 
overdensity is (1-H^(5l)^(1 + (5l)^ ~ (l-hff^i) (the local 
curved universe growth factor multiplied by a local-global 
conversion factor) . 

Combining these two effects ~ non-linear contraction 
and growth - we can write down the correlation function 
on small scales in the presence of a long wavelength mode 
Sl 



(1 + f <^L)e5o((l + y)r) 

e5o(r) + (-e5o(r) + -re^joM)'5L 



(2) 



68<5i Sl 



21 3 



where we denote linear quantities with subscript zeros 
and we have neglected other second order terms as well 
as other effects. 

As an aside, multiplying \>y Sl and averaging over re- 
alizations of the long wavelength modes we obtain an ex- 
pression for the squeezed bispectrum as an intermediate 
step: 



{51) 



(3) 



where S\ is the variance of the long wavelength mode. 

Does the presence of long wavelength modes shift the 
BAO feature to a different average scale? To find the 
shift predicted from our expression ([2]) , we must compare 
our expression with a shifted correlation function ^(r) = 

[ar) . Expanding for small shift values a — \ gives 



(4) 



Though in simulations the shift a — 1 is typically ob- 
tained from a least-squares fit, we can obtain an approx- 
imate estimate of the shift which best fits our expression 
^ by isolating the coefficient of ?'Cso('') after averaging 
over long wavelength modes. The lowest order non-zero 



coefficient is of order S\ ; 



68 , . 



63 



(5) 



We can see that in this simplified picture (of equation[2]), 
the BAO scale shift arises from the following effect: Be- 
cause perturbations in regions with long-wavelength over- 
densities grow larger, such regions contribute more to the 
correlation Junction, and so we tend to predominantly ob- 
serve overdense regions when measuring the BAO scale 
from the correlation Junction. In overdense regions, the 
BAO scale shrinks because such regions locally behave as 
positively curved closed universes, so that we measure a 
smaller scale than predicted by linear theory. If we mea- 
sure the correlation function using biased tracers, the 
effect should be larger, as we are even more sensitive to 
overdense regions and hence should see an even larger 
shift. 

Our derivation assumes that the BAO shift is due to 
the infiuence of long wavelength modes. Though this 
makes sense physically (short wavelength modes fluctu- 
ate rapidly over distances smaller than the BAO scale and 
hence should not lead to coherent gravitational shifts of 
it), we examine this assumption in the following section. 
Furthermore, our simple model neglects the influence of 
anisotropics (it assumes spherical symmetry) as well as 
other small corrections and neglects other possible second 
or higher order effects. A perturbation theory calculation 
including all these effects, leading to a slightly modified 
result, is shown in the following section. However, most 
of the shift can be described by the simple physical pic- 
ture just presented. 



III. PERTURBATION THEORY CALCULATION 
OF THE BAO SCALE SHIFT 

In this section we will use second order Eulerian per- 
turbation theory to derive a more rigorous expression for 
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the shift and compare it with our heuristically derived 
expression. Though other calculations using the machin- 
ery of renormalized perturbation theory and Lagrangian 
perturbation theory have been performed in |19l , |20| , a 
simple Eulerian approach is more pedagogical and still 
captures most of the physics (see also the similar discus- 
sion of the largest of the shift terms in [l9|). 

In standard perturbation theory we expand the density 
field in powers of the primordial perturbation Sq: 



S = 5(1) + ^(2) + ^(3) + 



(6) 



where S^"'^ contains i5o to the nth power. 

We will be interested in the power spectrum on the 
scale of the BAO feature ( "small scales" for our purposes) 
as well as the squeezed limit of the three point function. 
The small scale power spectrum has two contributions, 
(5(2)^(2)) and ((5(1) (5(3)); however, the 1-3 term can eas- 
ily be shown to be a simple rescaling of the linear power 
spectrum and hence does not contribute to the shift. Fur- 
thermore, as the three point function involves {dd^^^S^^")) 
to lowest order, we only need to know S up to second 
order to evaluate its leading order contribution. Hence 
to determine both the shift in the small scale power spec- 
trum and the squeezed Bispectrum we only need to study 
the structure of (5(^). 

In Fourier space (5(^-' is given by 



,5(2) (k) 

where 
i^2(k,q) - 



:^^2(k-q,q)^o(k-q)(5o(q) (7) 



17 1 k_q k 

21 2 kq q k' 



(8) 



We will initially work in real space for clarity, where 
the equivalent expression is: 

(5(x) = ^o(x)+d(x) . V5o(x) + ^5l{^) + ^i^,,(x)ify(x) - 
^ „ ' 2.1 _ 7 



shift 



where 



growth 



anisotropy 



(9) 



)3 q2 

^{'-^-h6o{ci)e'^-. (10) 
{2Try q'^ 6 



Here the shift term encodes the motion of the density per- 
turbations due to the gravitational potential sourced by 
other perturbations, the growth term describes the non- 
linear growth of perturbations due to the gravity from 
other perturbations, and the anisotropy term describes 
the anisotropic distortion of the perturbations. 



Hence we obtain to fourth order in 
,5(xi),5(x2) = i^V,5o(xi).d(xi),5o2(x2) + (1^2) 

+ yV5o(xi) • d(xi)if„(x2)if,,(x2) + (1 o 2) 
+ V(5o(xi)-d(xi)V(5o(x2)-d(x2) 

+ il^,52(xi)if„(x2)i^.,(x2) + (1 O 2) 



- ) K^,{yL2)K^j{y.2)K,j{:>^l)K^J{^l) 



+ ,5o(xi)(5(x2) + (1 o 2) + 



(11) 



where we have neglected the 1-3 term which does not 
contribute to the shift. 

We now take the average of this quantity so that we 
can derive the shift by identifying the coefficient of r^' (r) . 
To begin, we consider the first term from the equation 
above (the growth-shift term), at a radius of order the 
BAO scale rs'- 



(^V5o(xi)-d(xi)(52(x2) + (lo2)) 
21 



(12) 



^r^0{r) ■ J J2^~^^^^ ' ~ ^^^^ 



We can see that for modes with wavenumber q » r = 
xi — X2I ^ rg, the integrand above is highly oscillatory 
and averages to zero. This makes physical sense, because 
short wavelength modes cannot coherently shrink the cor- 
relation function on scales ~ rg much larger than their 
wavelength, as their gravitational effects average to zero 
over the BAO scale. Hence one can cut off the integrand 
for wavelengths above (/ > A, where A ~ l/j's, and ex- 
pand for qr < 1, so that the integral simply becomes an 
average over long wavelength modes (longer than r^): 



;i^V5o(xi)-d(xi)(52(x2) + (lo2)) 



-£'(r)A 
21 ^"^"^^ar 



68 
63 



dq 



<Pq Pojq) , 
(27r)3 g2 

Anq^Poiq) 



68 
63 



(Si) 



(13) 



where fi — ci ■ r. Though there is still a contribution to 
the integral from scales somewhat smaller than the cutoff 
A = I/t'b, this contribution does not scale as r and hence 
the term associated with it does not have the form of a 
shift of the correlation function ~ r(^'{r). Partially due 
to the approximate nature of the cutoff choice, it should 
be noted that our expressions for the shift can only be 
order of magnitude estimates. 
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All the other terms are similarly oscillatory (see Ap- 
pendix |Xl where they are all written in terms of Bessel 
functions) so that they can be similarly cut off. We thus 
obtain for each of the terms (ignoring a uniform shift 
independent of position as well as the 1-3 term); 



Cs(r) = (5s(xi)(5s(x2)) 



68 
63 
32 
315 



{51) 



(Si) 



15 



2 
7 

34 
21 



10 



2^50 (r) 



iso{r){6l) 



or succinctly 



,2438 



735 



131 
105 



10 



(15) 



Note that this expression can be obtained more directly 
by writing 5 = 8^ -\-5^ and keeping the long contributions 
to d, 5 and K but not to the derivative of 5 to give the 
expression 

5s{^) = ^5o(x+d^) + ^<5io'55o(x)+^i^|(x)i^,,^(x)+. • • , 

(16) 

and calculating the short mode correlation function 
^s(r) = {5s5s). 

As before, we identify the coefhcient of f^so('') in equa- 
tion (IT5|) as the shift of the real space correlation function 
a -I: 



105^ ^' 



(17) 



The 131/105 coefficient in front of the shift term 
is 131/105 = 68/63 -I- 32/315 -I- 1/15, where the first 
term arises from the cross between the shift and growth 
terms, the second term arises from the cross between the 
anisotropic and shift terms and the third term comes 
from the square of the shift term alone. The first term, 
which is what we obtained from our simple physical ar- 
gument, is the dominant term (the other two terms are 
less than 10% of its magnitude). 

In Fourier space our expression leads to: 



Ps{k) = Pso(fc) 



(k) - —kPUk) + —k^P''{k) 
735 ^ ^ 105 ^ 10 ' 



(Si) 
(18) 



where we have used 



-(fcP'(fc) -f 3P(fc)). (19) 



and 



dr 



rCir), 



d 
dka 



d 



okb 



12P(fc) + 8fcP'(fc) + k^P"{k). 



(20) 



The direct calculation in Fourier space (expanding such 
that one mode q has a much longer wavelength than the 
other, k » q) gives: 



Ps{k) 
(14) 



Psoik) 
d^q 



Psoik) + 
569 



d^q 



(2^)3-4[^2(k - q, q)]2Po(|k - q|)Po(g) 



735 



Psoik) 



47 
105 



10 



105^^^^ 



(22) 



{5LUr)) = [^isoir) + -resoirWl). 



(23) 



in agreement with our intuitive arguments. To go to 
Fourier space we transform f^sol'') give: 



^21 



Psoik) ~-kP's,ik)]{5l). 



(24) 



(51) 



(2.)3^°(^)^^^«(^) 



which is the same as our previous expression except for 
the fc^/3g^P(fc) term, which comes from the uniform 
shift we ignored in the real space calculation. This shift 
moves both points in the correlation function by the same 
amount so it leads to no effect. In the standard pertur- 
bation theory calculation this term cancels when P22 is 
combined with P13 (see Appendix VIII). 

The Fourier space shift ak describes how much the 
linear spectrum is rescaled as Psoik /ak)- Expanding in 
ak — 1 as before, we find we can determine it by identi- 
fying the coefficient of the shift term —kP'ggik) 



Note that the shift as calculated in Fourier space is dif- 
ferent from the shift in the real space correlation func- 
tion. This is due to the the second derivative terms mix- 
ing shift and non-shift terms upon Fourier transforming; 
physically, there is a change of BAO peak shape which 
looks like a shift in Fourier space. 

The squeezed limit of the three point function is 



which agrees with a direct Fourier space calculation of 
the squeezed Bispectrum. 
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IV. THE EFFECTS OF TRACER BIAS ON THE 
MEASURED BAO SCALE 



VI. THE EFFECT OF RECONSTRUCTION 
USING THE ZEL'DOVICH APPROXIMATION 



In order to compare our calculation with the semian- 
alytical results of [20|, we now include bias in the same 
way using the simplified expression 



Sh - bi6 + h}26^ + ■■■ 
The averaged two point function now gives 



(25) 



.2^r^2.2438 

,68 _ . 2 



f3f 



Previous measurements of the BAO scale have em- 
ployed reconstruction techniques to remove large scale ve- 
locity fields which broaden the BAO peak. These meth- 
ods [22] use the Zel'dovich approximation to attempt to 
restore the linear density field and hence undo the broad- 
ening and degrading of the BAO feature. In this section, 
we will briefly review the reconstruction procedure and 
examine to what extent reconstruction also reverses the 
non-linear shift of the peak. 

The Eulerian particle position x is related to the La- 



S,SQ{r) + _|^,.r^5Q(r) + ?5'o(''))grangian position r through the displacement vector q 



105 



10 



(26) 



x = q 



(32) 



In Fourier space this gives: 



As the small mass elements are related by pod^r = pd^x, 
we can relate the overdensity to the Jacobian and expand 
to first order: 



■735 



1 

To' 



26 2 
+ bMi^Ps,{k) - -kP's.m 



bjPsoik)}. (27) 



^ = 1 
Po 



d^r 
d^x 



We can hence identify the BAO scale shift in Fourier 
space, including the effects of bias, by finding the coeffi- 
cient of -bjkP'ik): 



47 , 

ttfe - 1 = 1 

105^ 



70 62 
47 bi 



(28) 



As 70/47 = 1.489, this expression agrees with the numer- 
ical calculation in [20| which gave at — 1 oc (1 -I- 1.5j^). 

The squeezed three point function can also be calcu- 
lated 



('5^^5.(0) 



ao 1 

{bl[^^S,{r) + -ri's,{r)]+2b,b2is,{r)){5l). 

(29) 



So that 



_ = J-i«l-V-g (33) 



-V • q (34) 



For the linear growing mode, q is curl-free so that one 
can define a potential q = —Vcfiz- In order to only undo 
large scale velocity flows, one smooths the density field 
(typically with a Gaussian of width ^ 20 Mpc) to obtain 
a new field 5'^ and then solves Poisson's equation for (f>z 
to obtain: 



d^q 1 
(27r)3^ 



<5^(q)e 



(35) 



Hence to undo the large scale velocity flows, one simply 
displaces the measured density field 5 by q = —^cfyz 
where 



V. BAO SCALE SHIFT: COMPARISON WITH 
SIMULATIONS 

Evaluating (5\) as the variance at the BAO 



cr2(i40Mpc) [D{z)/D{Q)Y gives 



scale 

Qm^l\D{z) / D{{))Y using WMAP cosmological parame- 



(30) 



ters [25|, so that from equation ([22 



ttfe - 1 w 0.23% [D{z)/D{Q)Y 



Considering that the numerical value depends on the 
choice of cutoff scale and window function, and that we 
have neglected J^a, this a gree s well with the value ob- 
tained from simulations in llTl of 



ak{z) - 1 = (0.30 ± 0.02)% [D{z)/D{0)f 



(31) 



d^q i<\ 

(27r)3 ^ 



<5G(q)e 



zq-x 



Comparing this with the shift 

d^q zq 
(2^^ 



^Lo(q)e* 



(36) 



(37) 



in equation (17), which relates the linear and non-linear 
fields as 



5s(x) = <5so(x + d^) + -5lo<55o(x) + -A;^-(x)i^,^(x), 

(38) 

we can see that q is (very nearly) identical to — and 
hence reconstruction (i.e. displacing the field by q) will 
undo the non-linear shift we found previously, as long 
as it includes in 8^ all the long wavelength modes re- 
sponsible for the shift. As the broadening is due to a 
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wider range of modes than the shift (including shorter 
wavelengths), if a reconstruction procedure reverses the 
peak broadening, it automatically also reverses the low- 
est order shift. Our conclusion that the BAO shift should 
be greatly reduced after reconstruction agrees with sim- 
ulation results from [23|, who found that the shifts were 
consistent with zero after reconstruction was applied. It 
also agrees with perturbation theory calculations in . 

VII. CONCLUSIONS 

We have investigated the non-linear shift in the BAO 
scale found in simulations and have tried to develop an 
understanding of why the scale becomes smaller due to 
non-linearities. Using both intuitive arguments and per- 
turbation theory, we found that the main reason for the 
shift in the BAO scale is as follows: As perturbations 
in regions with long-wavelength overdensities are larger 
due to non-linear growth and bias, these regions con- 
tribute more to the correlation function so that a mea- 
surement of the BAO scale is more sensitive to overdense 
regions; in overdense regions, the BAO scale shrinks be- 
cause such regions locally behave as closed universes, so 
that we measure a smaller scale than predicted by linear 
theory. 

We have expanded and verified these heuristic argu- 
ments using perturbation theory, and derived expressions 
for the real and Fourier space shifts: 

a-l-^(^i); (39) 
where the average is taken over modes with wavelengths 



longer than the BAO scale (and can be approximated as 
(61) ~ a^iUOMpc)[D{z)/D{0)f). We note that, sur- 
prisingly, the values of the shift found if fit for in Fourier 
space should differ from the shift measured in real space. 
Using a simple model for the shift of biased tracers, we 
obtained an analytic expression in agreement with |20| . 
For dark matter our expression for the shift evaluates to 
Q!fc — 1 ^ 0.23% [D{z)/ D{0)]^ , which is in agreement with 
simulations. We explained in our physical picture why 
reconstruction of the BAO peak should undo the shifts, 
as also found in simulations. The sub-percent level of 
the BAO scale shift and the ease with which it can be 
undone using reconstruction confirm that it will not be 
problematic for surveys in the foreseeable future. 
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VIII. APPENDIX: CALCULATION OF THE 
CORRELATION FUNCTION TO THIRD ORDER 



In Fourier space we will write S^^^ as: 

d^qi (Pq2 (fqs 



^ = y (2.)3(2.)3(2.)3 ^^(q^'q^'q^)^ 



(27r)3 (27r)3 (27r)3 

<5o(qi)<5o(q2)<5o(q2)(27r)35^(k - (qi + q2 + qa)), (40) 



where we have defined 



with 



^3(qi,q2,q3) = -P5(qi,q2,q3) + -P3(q2,qi,q3) + -p3(q3,q2,qi) 

^3(qi,q2,q3) = ■r-{2/3(qi,q2 +q3)G'2(q2 + qs) + 7a(qi,q2 +q3)-F2(q2 +q3) 
54 



a(q, k) 



+(2/3(q2 +q3,qi) + 7a(,q2 +q3,qi))G'2(q2 +q3))} 
(q + k) • q 



g2 



i-i \ 17 Ik-q ,fc q. 2r,k-q,2 1, 
^ /, N 13 Ik-q ,fc q. 4r,k-q,, 1, 
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These are the standard perturbation theory kernels (see for example [2g , although we have rearranged some of the 
terms and thus we have introduced ^3. In our notation i^2j G2 and are the standard symmetrized kernels. 
One can use these kernels to compute the 1 — 3 contribution to the power spectrum in the limit q << k, 

/^3„ f j3„ 1 o-|l2 2 



116^ /M/x2\ f d'^l D ^ 



Psoiml) - / ^Po(q) Po(fc) (43) 



Note that as discussed before the piece proportional to 1/q^ cancels in the sum P22 + -Pis to give: 

P ^ P22 + Pi3 - i^Psoik) - ^kP'soik) + ±k^P'^,{k)}{Sl) (44) 

We can also use the above expression for (J^'^^ to go directly to real space. We are interested in the effect of a long 
mode on the short fluctuations and want to use iJ^'^-' to obtain the piece that is linear in the short modes and quadratic 
in the long. We obtain: 

(5(3)(x) = idf4a,c',5^(x)+A,a,<5^(x)+P(5^(x) 

+ a,Kf^{^) + D,,km%{^) + E,,,kiKl^i{^). (45) 
We have introduced A^, B, Cij, and Eij ^i which are all quadratic in the long mode and are given by: 

^ - 657^^ ) +189^^^*'^ + 27^»^-^^-'- + 18^ +54^ 



189 9 ''^ 9 

21' 



E^,M = ^S^Kt^Kti. (46) 



For convenience we have introduced: 



KatM^) = I - y^^)^^ - -^ScMci)e'^-. (47) 



The second order fluctuations are given by: 



'^('''^(k) = y'^F2(k-q,q)5o^(k-q)5„^(q) 
^^'^"-(k) - / ^G2(k-q,q)Jo^(k-q)<5o^(q) 

q 

Kb (q) = (— - n^ab)S'o (q). (48) 
q o 



In real space: 



5(^)^(x) = 4d,6^+^2^6^r + ^Kt,K^^ 

0(2)^(x) = 4d,S^ + ^{6^r + ^Kt.Kr (49) 

If we are interested in using these expression to compute the monopole of the short scale correlation function, 
then the terms proportional to K^^. will not contribute so we can ignore the Cij and Dijk terms. The terms 
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l/2dfdj'didj6^{yi) + Aidi6^{x) are responsible for a shift of the small scales that cancels in the correlation func- 
tion. Thus only the the terms proportional to B and Eij,ki contribute. Note that Kf^ cannot be neglected, and 
when computing averages over angles should be replaced by: 



Kfj^ki ^ [-^S,j6ki + —{5,k6ji + 5a6kj)]6^. (50) 



1 

45 Y5 

Thus for the purpose of computing the angled averaged correlation function we can use 



When calculating the correlation function, the contribution {S^^^S + SS'^'^^) leeds to a term we can call ^fg given by: 



„ ,682 208 116 2. 2^ , , 116 , , 

2 X 1 WHoir) = crHoir) (52) 

^567 189 945 3^ ^^"^ ' 315 ' ^ ' 



We can use equation p9|) and (|5ip to obtain: 



which when calculating the correlation function gives: 



= 2xgaiC5o(r) = ^ai^5o(0 (54) 



IX. APPENDIX: SIMPLE 5^ 

The full expression for 5'^ is easy to obtain in real space: 



(55) 



where 



kC{^) ^ /(03(^-^^'=™)^*'Hq)e--. (56) 

X. APPENDIX: 2-2 CORRELATION FUNCTION IN REAL SPACE 

The second order expression for the density field can be used to directly obtain the 2 — 2 piece of the correlation. 
We obtain: 

C(0|2-2 = ^Coo 

68 2 , 136 68 2 

+ ^Co,2C-l,l + Y^C0,2C0,0 + ^Co.O 

3 1 1 1 32 I2 

+ 2''2,2C-2,2 + 2C2,2C-2,0 + 2C2,oC-2,2 + 2C2,oC-2,0 + ^^0,2 + C2,oCo,0 + 2Co,0 

4 12 12 20 

— gCiaC-l,! ^Cl,3C-l,l ;^C-l,3Cl,l y-Cl,3C-l,3 

+ ^2+^0^^ + 5 2 + 20 ,6 11 2 

+ y^.0,2 + ^S0,4t,0,03 + yt,04 + Y^<.0, 21,0,0 + ^t,0,4<=.0,0 + ^(,0,0- \^ ' ) 
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The different lines correspond to tlie different terms in the calculation, growth-growth, growth-shift, growth- anisotropy, 
shift-shift, shift-anisotropy and anisotropy-anisotropy. 
We have defined 



Using this notation 



Cm,n — / 


dlnfcA2(fc)fc"[— 
okr 


Co,o(r) - 


Cir) 


Ci.iW = 


Cir) 


C2,2(r) = 


Cir) 


C2,o(r) = 


-v^CM = -C"M 



(58) 



(59) 



r 

We can now take the limit when one of the momenta is much smaller than 1/r which gives: 

dC(r)|2-2 1156^ 



dal 441 



■Cso 



68 , 

+ 

- 3?^^^° + T0"^^°+15"^^°+3^^° 

+ 315''^^° 
64 

Cso, (60) 



2205 



which adds to 



2 
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